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Abstract 

The entropy of a spherically symmetric distribution of matter in self- 
equilibrium is calculated. When gravitational effects are neglected, the en- 
tropy of the system is proportional to its volume. As effects due to gravita- 
tional self-interactions become more important, the entropy acquires a cor- 
rection term and is no longer purely volume scaling. In the limit that the 
boundary of the system approaches its event horizon, the total entropy of the 
system is proportional to its area. The scaling laws of the system's thermo- 
dynamical quantities are identical to those of a black hole, even though the 
system does not possess an event horizon. 
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The entropy of ordinary matter is usually proportional to its volume and depends cru- 
cially on its composition. The entropy of a black hole [|T| [|2[ is therefore very mysterious, in 
part because it is proportional to the black hole's area A. In order to explain this peculiar 
behavior, researchers have invoked holography, string theory, entanglement entropy , and 
brick-wall models 

However, the volume scaling properties for the entropy of ordinary matter depend on the 
fact that interactions are short range in comparison to the size of the system, or that the 
interactions are screened. For long-range interactions such as electromagnetism, the entropy 
still scales as the system's volume as long as the net charge of the system is neutral [^, 
but little is known about the thermodynamic properties in other cases where long range 
interactions are present. 

In this paper, we will examine a system where long range gravitational interactions are 
present. Although black hole entropy appears to possess very peculiar properties, we will 
see that many of these features are present in ordinary self-gravitating systems which do 
not possess an event horizon. In particular, we will study spherically symmetric matter 
distributions which are in self-equilibrium. In the absence of gravity, the entropy of the 
distribution is assumed to be proportional to its volume. However, when the system's own 
weak gravitational field is taken into account, the entropy acquires a correction term, spoiling 
its volume-scaling properties. 

In the strong field limit, the results are even more surprising. In the limit that the 
system is about to form a black hole, its entropy is proportional to its area for a large 
class of equations of state. The scaling laws of the system's temperature and energy are 
also identical to those of a black hole despite the fact that no horizon is present. This 
suggests that many of the peculiar properties of black hole thermodynamics are the result 
of gravitational self-interactions, rather than the horizon. 

In this work, the system is modeled as a dense series of thin shells of arbitrary compo- 
sition. This allows the system to be compressed to a size close to its own event horizon [Q. 
The work in this paper is related to reference where a single shell is considered. In that 
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case, the entropy is assumed to always be proportional to the system's area, and when the 
shell is brought close to its event horizon, the constant of proportionality is found to be 1/4. 

We consider n densely packed shells, where the radius of the i'th shell is given by r,, 
and its thermodynamical quantities, as measured by local observers on each shell are given 
by Ei (mass-energy), Tj (temperature). Si (entropy). Pi (tangential pressure), iVj (particle 
number), and /Xj (chemical potential). The entire configuration is assumed to be in equilib- 
rium with itself, and therefore, it supports itself in its own gravitational field. If the material 
has a rest mass, then this corresponds to a constant term in the chemical potential. One 
could also consider an arbitrary number of particle species, however the effect of this on our 
calculation is trivial. 

We assume that in the absence of gravity, the system has no unscreened long-range 
interactions and is therefore extensive. In other words, the entropy, particle number and 
energy scale as the size of the system when the intensive variables (temperature, pressure and 
chemical potential) are held fixed. In most applications of thermodynamics this is almost 
always assumed to be the case (at least implicitly), for without it, taking the thermodynamic 
limit becomes problematic. 

We will first calculate the local thermodynamical quantities in terms of observables at 
infinity. Spherical symmetry dictates that the metric outside the i'th shell is 

ds^ = iQ^j^dx^dx" 

= -Cik^{r)df + k-\r)dr^ + r^^dO^ + sin^ Odcj)^) (1) 

where throughout this paper we use units G = c = h = kB = l. In order for the time 
coordinate t to be continuous throughout the space time, the Cj need to satisfy the condition 

Ciki{ri+i) = d+ik'^^^in+i) (2) 

(an overall constant is not physically relevant, and just scales the time coordinate). If iT^t 
is the intrinsic stress-energy tensor on the i'th shell, then its components are related to the 
junction conditions of the extrinsic curvature and intrinsic metric [§] 
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SniTab = i[gabK - Kab] (3) 

where a = 0..2 and the notation J ] indicates the jump in the bracketed quantity evaluated 
at Tj e.g. = «^(^«) ~i-ik{ri). Since the shells are static, the local energy Ei of the i'th 

shell is —AnrfiTl, while the tangential pressure Pi is jT^ = jTg. The jump condition thus 
give 

= -^[k{r)]n (4) 
k{r 



SttR 



k'(r) + 



(5) 



If we take the material to have zero net charge (as indeed we must, for the system to be 
locally extensive), and have no angular momentum, then the ki{r) are of the Schwarzschild 
form 



fc.(r) = ^l-^. (6) 

Using (|^) we can then find ig^^ by calculating the constants given in 

^ _ (r^+i - 2mi+i)(ri+2 - 2m^+2) • • • (^n - 2m„) 
(ri+i - 2mi){ri+2 - Sm^+i) ■ ■ ■ (r„ - 2m„_i) 

(with c„ = 1). 

The energy, and pressure of each shell are then 



Ei = y/niyri - 2mi_i - Vn- 2mi) (8) 
p _ mi_i - Ti _ rrii - Vj 

The latter equation can be thought of as the condition for mechanical equilibrium, and 
ensures that the pressure balances the effect of gravity so that the configuration is at least 
momentarily static. The conditions for thermodynamic and chemical equilibrium are the 
Tolman relation PI 
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-J- r\ 



(10) 



-i9tt[ri) 

and the equality of the red-shifted chemical potentials 

l^iV i9tt — f^jy/j9u ■ 



(11) 



The latter ensures that there will be no net particle flow (one can imagine that particles 
can hop from shell to shell via wires, or simply by jumping). The Tolman relation fixes the 
temperature of the shells with respect to the temperature at infinity Too which we assume 
to be finite (to be rigorous, one could imagine a reservoir at infinity which is attached via a 
heat-conducting wire to all the shells). 

These thermodynamical variables can be related to each other by the Gibbs-Duhem 
relation |jlO| 



Ei = TiSi - AiPi + f^iNi 



(12) 



where Ai is the area of each shell. This equation holds for systems which are locally extensive. 
To derive this, we hold the intensive variables {Pi, Ti, /ij) fixed, and note that if we scale 
the size of a shell by A, the entropy and energy will scale by the same amount. I.e. 



Ei{XSi,XAi) = XEi{Si,Ai] 



(13) 



We can then differentiate the above expression by A to obtain the Gibbs-Dehum relation 
(which is an Euler relation of homogeneity one). In the presence of gravity, the equivalence 
principle ensures that the Gibbs-Duhem relation is valid as it applies to thermodynamical 
variables as measured in the proper rest frame of the shell. These variables are scalars, and 
are therefore frame- independent. 



Substituting the expressions for the local variables into (|T2D one can write the entropy 
of each shell as 



S,. 



1 
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{3mi - Ti) - (3mj_i - r^) 



ki-i{ri) 



'c, - 



T 



'1 



2mr 



r„ 



(14) 
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and the total entropy of the system is just the sum 

S = T.S^ . (15) 
i 

When gravitational effects are negligible, i.e. the configuration of shells is large in com- 
parison to its Schwarzschild radius, Si can be expanded to first order in rrii/ri. It is instruc- 
tive to look at some of the terms in equation ( |12|) individually. The local energy can be 
expanded to give 

i?, ^ - m,_i + ~ ""^-^ . (16) 

rrii — mj_i can be interpreted as the additional energy (as measured at infinity) required to 
add the i'th shell to the configuration. Since to first order, mj is just the sum of all Ej, 
j < i, v^e see that we have recovered the usual Newtonian expression 

rrii - mi_i = Ei + (pi (17) 

where (pi is the Newtonian gravitational potential 

0. = _E. ^^ + ^^ + -^-^ _ f . (18) 
ri 2ri 

To first order in rrii/ri, the pressure term also gives the usual Newtonian expression 
The chemical potential term 



fxjl 1-— 20 

V rn r„ 

contains a correction which is equivalent to the Newtonian chemical potential due to an 
external gravitational field, although in this case, the gravitational field couples not only to 
the mass of each particle, but to the entire chemical potential fin- 
One can also expand the expression for the temperature, and after substituting these 
terms in (^4|) and taking the sum over i, we arrive at the total entropy. To zeroth order 

6 



Sc^T-\mn~^^nN) (21) 

where is the total number of particles. We see that when we ignore the first order effects 
due to gravity, the entropy, when expressed in terms of observables at infinity, is given as an 
Euler relation of homogeneity one. This demonstrates that the system is purely extensive. 
When the intensive variables (Too and are held fixed, the entropy is proportional to the 
total energy and total number of particles. If all the shells have the same energy density 
and number density, the entropy would scale as the volume of the system, where the volume 
is understood to be the sum of all the Ai. 

However, if we include the first order terms, we find 

(b Tfi 

S - T-\m^ - /i„iV + ^ - /i„iV— ) (22) 
where (p is the total gravitational energy required to construct the system 

= ^0, . (23) 

i 

The entropy includes correction terms, due to the gravitational self-interactions, and is no 
longer purely extensive. 

We now go to the extreme limit when the outer shell approaches its own Schwarzschild 
radius (r„ ^rrin)- For the purposes of this paper, we will assume that our system can 
approach its Schwarzschild radius, although more physically reasonable models may not 
be able to generate such intense pressure or exist at such high temperatures. Inspecting 
equation (0) we see that all the q except c„ are zero because of the r^r — 2mjv term. As a 
result, from equation (11^) we see that all the Si except Sn are zero, although for now, we 
will retain the chemical potential terms. Substituting this result in (^) we can easily sum 
over i and calculate the total entropy 



^ = ^-^Jl-^ . (24) 



2T TV r 



Pretorius, VoUick and Israel consider specific and general equations of state where the 
chemical potential term in equation ( P^ vanishes entirely However, we can also demand 
the less restrictive condition that the chemical potential not be overly divergent i.e. 



/ 2m 

lim /i„A^Jl -<^mn . (25) 

As a result, the only term which contributes to equation (|15]) is the entropy of the final shell 

S = ^ . (26) 

Since the total entropy is equal to the entropy of the final shell, we see that the entropy is 
proportional to the area of the system. That this is the case, can also be seen by noting that 



equation (^61) is an Euler relation of homogeneity two. This is the precise relation which 
governs the entropy and temperature of black holes. Applying the first law drrin = T^dS to 
(^), one finds that (in terms of a constant 7) 

S = 'yml 

^ An (27) 



167r 
and 

Too = (27m„)-' (28) 

where as before, m„ is the total ADM mass as measured at infinity, and An is the area of 
the final shell (and therefore, of the entire system). 

These relations, up to the constant 7, are identical to the ones which govern Schwarzschild 
black holes in 3 + 1 dimensions. Here however, the entropy is not a property of a horizon, 
but is the logarithm of the number of micro-states of our system. The temperature, which is 
usually considered to be independent of the size of the system, is now inversely proportional 
to the mass of the system. This result is even more remarkable, since it has not been derived 
by considering quantum fields on the space-time. These later two equations are the central 
result of this paper. 

It is also interesting to note that if, at the instant before the system forms a black hole, 
it has the Hartle-Hawking temperature 
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Thh = 2m„/A 



-n 



(29) 



then equations (^) and ( pTf ) imply that its entropy before black hole formation is 



S{THH,mn) = An/ 4: 



(30) 



This is identical to the entropy of a black hole of equivalent mass. 

Indeed, Pretorius, Vollick, and Israel argue that equation (^) is the appropriate temper- 
ature for a shell as r„ 2m„. They consider quantum fields propagating on the space-time, 
and demand that that the shell be in equilibrium with the local acceleration temperature 



Ordinarily, the appropriate vacuum state for a spherically symmetric distribution with- 
out a horizon is the Boulware vacuum (which has a temperature of zero at infinity). However, 
when the boundary of our system lies close to its Schwarzschild radius, one cannot consis- 
tently construct the Boulware vacuum, since the stress-energy of the quantum fields will 
be infinite (and negative), and will depend on the number of fields. These trans-plankian 
frequencies will result in an enormous and incalculable back-reaction. It is therefore unclear 
what the appropriate vacuum state and vacuum temperature are. Pretorius, Vollick, and 
Israel demand that the quantum field be in the Hartle-Hawking state with a temperature 
at infinity given by equation (^), however, this requires adding layers of insulation to the 
shells, and positive energy to the field to compensate for the back-reaction. 

The manner in which this is done, and the resulting temperature of the shells is not 
unique. For example, for a single shell just outside its Schwarzschild horizon, consistently 
constructing the Hartle-Hawking state requires adding a thin layer of insulation to the shell. 
The insulation prevents the fields inside the shell from having a diverging temperature, which 
would result in an infinite back-reaction. If the insulation is added to the outside of the shell, 
the temperature of the shell can be zero, since it is shielded from the hot field just outside 
the shell. If the insulation is added to the inside of the shell, the local shell temperature will 
match the local temperature of the field and will therefore diverge as Vn 2mn- 



0- 
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The resulting scaling laws we have derived apply when the system is in static equilibrium 
just outside its Schwarzschild radius. By extension, they also apply at the point right before 
the system forms a black hole, if the collapse is quasistatic (i.e. the system is at equilibrium 
at all points during collapse). Any system which forms a black hole through a collapse 
off of equilibrium, will have an entropy bounded from above by equation (P7D, since the 
equilibrium configuration maximizes the entropy. Furthermore, for a system which forms a 
black hole, the generalized second law |l|] implies that 7A„/(167r) is bounded from above by 
A„/4, since otherwise it would not be entropically favorable for the system to form a black 
hole. 

We have seen that the usual volume scaling properties of a system no longer hold when 
gravitational effects are taken into account. In the limit that our system approaches its 
Schwarzschild radius, its entropy scales as the area of the system. The temperature is also 
no longer an intensive variable, and scales inversely with the system's total energy. In other 
words, the system's thermodynamical quantities scale in the same manner as those of a 
black hole. These rather surprising results suggestions that at least to an extent, some of 
the intriguing aspects of black hole entropy can also be found in self-gravitating systems 
which do not possess a horizon. Once long range interactions are taken into account, the 
entropy of a system can behave in unfamiliar ways, and lead to scaling behavior which is no 
longer extensive. 

The ability of high-energy theories such as loop quantum gravity and string theory to 
make predictions about the black hole entropy has generated much interest, and provides 
crucial insights into the entropy of black holes. The results in this work suggest that a wide 
class of systems will also result in the correct entropy scaling laws for black holes. It is hoped 
that by studying the entropy of these self-gravitating systems, we can learn more about black 
hole entropy. It is also hoped that studying the thermodynamics of self-gravitating systems 
will lead to a better understanding of the thermodynamics of more general non-extensive 
systems. 
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